Biomechanics and Modeling in Mechanobiology

The mechanical behavior of right ventricular (RV) myocardium is governed by its anisotropic microstructure, yet constitutive models that account for fiber dispersion and enable reliable parameter identification remain limited. In this study, we propose a physics-embedded constitutive neural network framework for automated discovery of strain energy functions and microstructural parameters from experimental data. The model is formulated within an incompressible, orthotropic hyperelastic setting using invariant-based representations. Fiber, sheet, and normal directions are incorporated through a rotated structural basis, and dispersion effects are modeled using a generalized structure tensor approach. The framework is trained on multi-axial mechanical data from ovine RV myocardium, including uniaxial tension-compression and simple shear tests. We investigate two training scenarios: (i) full datasets containing both tensile and compressive regimes and (ii) datasets restricted to tensile loading. In both cases, the model accurately reproduces the measured stress-strain responses and identifies sparse, interpretable constitutive models which involve isotropic, anisotropic, and coupling invariants. However, the identifiability of microstructural parameters strongly depends on the available loading conditions. While tensile-only data yield higher predictive accuracy, they result in non-unique or biased estimates of fiber dispersion. In contrast, inclusion of compressive data enables consistent identification of dispersion parameters by separating fiber and matrix contributions. These results highlight the importance of multi-axial loading data for robust parameter identification and demonstrate the capability of constitutive neural network-based approaches for data-driven modeling of anisotropic soft tissues.

The structural integrity of spinal cord tissue and the transmission of mechanical stimuli across the different levels of tissue microarchitecture and varying spatial scales of mechanical loading challenge experimental and computational efforts to accurately model, simulate and interpret tissue mechanics, leading to conflicting findings in existing literature. Here, we demonstrate that the bead size used in spherical indentation tests significantly affects the stiffness ratio of spinal cord gray to white matter, a dependence which we only observe on the transverse plane and not the coronal plane of the tissue. Our study reveals a shift in stiffness ratio such that for smaller spherical indenters gray matter is stiffer than white matter, while for larger indenters, white matter is stiffer than gray matter. The mean relative change from the 100 {micro}m bead to the 500 {micro}m bead differed between anatomical planes, with transverse sections showing a decrease in gray matter (-13.3%) and an increase in white matter stiffness (+26.9%), accompanied by a reduction in the gray-to-white matter stiffness ratio from 1.07 to 0.76, whereas coronal sections exhibited increases in both gray (+21.0%) and white matter (+33.8%), along with a change in the ratio from 0.99 to 1.14. These findings contribute to explaining previously contradictory results in the literature and underscore the relevance of spatial scales in mechanical characterization studies.

Finite element (FE) head models are valuable tools for investigating brain injury mechanics, with their reliability critically dependent on accurate material modelling. White matter (WM) is often considered mechanically anisotropic due to its aligned axonal fiber architecture and is commonly represented using fiber-reinforced hyperelastic formulations such as the Gasser-Ogden-Holzapfel (GOH) model. A fundamental assumption of the GOH model is that fibers contribute only in tension and not in compression, requiring the use of tension-compression switches. However, inconsistencies were noted in the formulation of tension-compression switches with the influence on computational biomechanics unknown. To address this knowledge gap, three commonly used switching schemes - differing in both the switching parameter and the treatment of compressed fibers - were theoretically elaborated and numerical implementation within the GOH framework to simulate the mechanical anisotropy of WM in impact simulations. Results from the case-based and group-level analyses demonstrated that both the switching parameter and the treatment of compressed fibers affected WM deformation. Significant cross-scheme strain differences were noted in the first principal strain at the element level and fiber strain at the fiber level. These findings highlighted the mechanical role of tension-compression switch in the GOH-based brain modelling and advocated the adoption of fiber stretch itself as the switching parameter to discriminate the tensile and compressive fibers. The current study provides important guidance for the anisotropic constitutive models in brain tissue and calls for direct verification of the tension-compression switch hypothesis in axonal fibers.

Functional spinal cord repair in zebrafish is governed by regeneration-favorable biochemical and mechanical cues within the lesion microenvironment. Alterations in extracellular matrix composition and stiffness are closely associated with axon regeneration. However, experimentally dissecting the interplay between mechanical signals and axonal regrowth in vivo remains technically challenging. Here, we present an agent-based modeling framework to simulate stiffness-mediated axonal growth trajectories across the lesion. We use this model to explore potential mechanisms underlying the characteristic growth patterns observed during zebrafish spinal cord regeneration. Computational predictions were qualitatively compared with confocal imaging data obtained from larval zebrafish. These phenomenological comparisons revealed a close agreement between simulated and experimentally observed axon growth, indicating that experimentally observed patterns could be governed by transient changes in the stiffness profile of the spinal cord and lesion microenvironment. Hence, our computational framework provides an in silico platform for investigating the role of mechanical cues in axon regeneration in the injured spinal cord.

The aorta, normally resilient to hemodynamic stresses, becomes vulnerable to structural failure due to diverse conditions that weaken the wall. We injected fluid into excised specimens of human ascending aorta with pressure monitoring to quantify the impact of clinical and histological factors on mural damage. Two modes of medial injury emerged with distinct pressure tracings. Extravasation was characterized by diffuse infiltration of fluid with widespread damage of smooth muscle cells and collagen fibers but limited separation of elastic lamellae. By contrast, delamination was characterized by marked separation of elastic lamellae along a single plane with damage to cells and fibrillar matrix restricted to adjacent laminae. Aging, aortic dilatation, and family history associated with lower pressures causing delamination, whereas a diagnosis of hypertension associated with higher pressures suggesting resilience to dissection. Collagen fraction adjacent to delamination correlated with higher pressures as did decreased smooth muscle cell density and increased glycosaminoglycan fraction, although several clinical and histological variables were interrelated. Protein cross-linking strengthened and enzymatic digestion of collagen weakened the aortic wall, while acute cell lysis with detergent had no effect. We conclude that increased functional medial collagen has an adaptive protective role in aortic remodeling rather than signifying medial degeneration.

Clavicle fractures often exhibit markedly different clinical outcomes: some patients recover acceptable function despite shortening or displacement, whereas others with apparently similar deformity develop persistent pain, functional loss, or poor healing. To explain this distinction, we propose a minimal nonlinear mechanical model for prognostic analysis of clavicle fractures. The model describes the interaction between fracture-related shortening and compensatory shoulder-girdle posture through a reduced equilibrium equation incorporating stiffness, geometric nonlinearity, and shortening-posture coupling. Within this framework, we analyze equilibrium branches, local stability, and the emergence of critical thresholds. We show that post-fracture destabilization can be interpreted as a fold bifurcation, while more complex parameter dependence gives rise to cusp-type structures and multistability. These bifurcation mechanisms provide a mathematical explanation for sudden deterioration after injury or treatment, as well as for strong inter-individual variability. We further introduce an optimization principle based on a utility functional to guide treatment planning. The analysis predicts that the optimal safe correction should lie strictly below the bifurcation threshold, thereby generating a natural safety margin. Although the model is simplified and has not yet been calibrated against patient data, it nevertheless provides a theoretical framework for understanding why fracture prognosis may deteriorate abruptly near critical mechanical conditions and offers a dynamical-systems interpretation of empirical treatment thresholds used in clinical practice.

Thoracic Endovascular Aortic Repair (TEVAR) is a minimally invasive procedure for the treatment of thoracic aortic pathologies, such as Thoracic Aortic Aneurysm (TAA). Computational simulations can provide valuable insights into TEVAR outcomes and complications prior to surgery, making them a useful tool in the procedural planning. In this work, Fluid-Structure Interaction (FSI) computational simulations are carried out in ten pre-TEVAR patient-specific TAA cases, for which post-TEVAR outcomes are known, to quantify the hemodynamic drag forces acting on the aortic wall. Based on these results, this study proposes a new risk factor R to predict the occurrence of type I and III endoleaks. The patient cohort is divided in a calibration set, used to associate specific R values with three different risk levels, and a validation set, to test the risk factor efficacy. Based on the risk factor values obtained for the calibration set, R[&le;] 0.33 is associated with low risk of endoleak formation, 0.33 < R[&le;] 0.67 with moderate risk, and R > 0.67 with high risk. Once it is applied to the validation set,the risk factor is able to predict the formation of a type Ia endoleak. The risk factor proposed in this work is capable of identifying all the endoleak cases analysed, as well as conditions known to increase the risk of TEVAR complications. This study represents a preliminary attempt to determine whether pre-TEVAR hemodynamics can effectively predict post-TEVAR complications and thereby aid clinicians in the pre-operative planning.

Thin structures such as heart valves and aortic dissection flaps interact dynamically with blood flow in human vessels. Their flexibility and capacity for large deformations generate complex, highly transient hemodynamic patterns over the cardiac cycle. Accurately resolving these interactions remains challenging for conventional boundary-fitted fluid-structure interaction approaches. We present an immersed boundary method for simulating thin structures in incompressible flow on unstructured grids. The method couples a stabilized finite element fluid solver with a nonlinear, rotation-free shell formulation through a direct forcing immersed boundary approach. The framework supports both weak (explicit) and strong (implicit) time-coupling strategies, enabling stable simulations over a wide range of solid-to-fluid density ratios. Hydrodynamic forces acting on thin structures are computed from fluid solutions sampled on both sides of the structure, allowing accurate force reconstruction for zero-thickness shells. To our knowledge, this is the first immersed boundary formulation that couples an unstructured finite element fluid solver with a two-dimensional, rotation-free shell model to simulate interactions between thin structures and incompressible flow. Fluid-structure coupling is achieved using predefined finite element shape functions, which provide consistent projection between Eulerian and Lagrangian fields without additional interpolation procedures. The framework is validated using three-dimensional benchmark problems involving thin structures. Then, valve-like model is used to compare strong and weak coupling strategies. Finally, the method is applied to an idealized type-B aortic dissection model. The proposed approach is implemented within the open-source software CRIMSON, a finite element platform for cardiovascular simulation.

Cerebral blood flow (CBF) control is essential for normal brain function and is disrupted in pathological conditions. Arterial diameters are tightly regulated to provide on demand increases in blood flow in regions of neuronal activity. Pericytes (PCs) exhibit robust myogenic tone and may also respond to neuronal activity to fine-tune local resistance and blood flow. Thus, mural control of microcirculatory resistance may extend beyond arteries and arterioles. Yet, PCs electrophysiology and contractility have not been thoroughly characterized, and this prohibits an integrated view of brain blood flow control. In this study, we develop a detailed mathematical model of mural cell electrophysiology, Ca2+ dynamics and biomechanics. The model is informed by electrophysiological data in smooth muscle cells (SMCs) or PCs and predictions are compared against pressure-induced responses in isolated arterioles and capillaries, respectively. Simulations recapitulate myogenic constrictions and examine differences in contractile dynamics as we move from arterioles to proximal and distal capillaries. In arteriole-to-capillary transitional (ACT) zone PCs, increased mechanosensitivity, more Ca2+ influx through non-selective cation (NSC) channels and/or a higher sensitivity of the contractile apparatus to Ca2+ can compensate for reduced L-type voltage-operated (VOCC) Ca2+ influx and allow for robust constrictions at the lower operating pressures of capillaries relative to the arterioles. A significant Ca2+ influx through NSC relative to VOCC, however, can decouple the PCs contractile apparatus from electrical signaling. Vasoactivity to chemomechanical stimuli along the arteriole to capillary axis is progressively driven by VOCC-independent Ca2+ influx and Ca2+ sensitization with slow kinetics. The proposed cell model can form the basis for detailed multiscale and multicellular models that will examine physiological function at a single vessel or vascular network levels and investigate CBF control in health and in disease. Key pointsO_LIA mural cell model of electrophysiology, calcium (Ca2+) dynamics and biomechanics is informed by data and adapted for modeling cerebral arteriole smooth muscle cells and capillary pericytes. C_LIO_LIIon channel activities are characterized by patch-clamp electrophysiology in isolated cerebral smooth muscle cell and pericytes, and capillary and arteriole electromechanical responses to transmural pressure changes are assessed using novel ex vivo preparations. C_LIO_LIMyogenic constrictions in arterioles can be reproduced by pressure-induced non-selective cation channel (NSC) activation that depolarizes the cell, opens L-type Ca2+ channels (VOCCs) and increases Ca2+ influx. C_LIO_LIRobust myogenic constrictions in arteriole-to-capillary transition (ACT) zone pericytes may reflect significant Ca2+ influx through NSC, increased mechanosensitivity, or higher sensitivity of the contractile apparatus to Ca2+, potentially compensating for reduced VOCC density relative to arteriolar smooth muscle. C_LIO_LIA significant contribution of NSC relative to VOCC in Ca2+ influx, can decouple the contractile apparatus from electrical signaling. C_LIO_LIThe model shows how gradients in ionic activities, mechanosensitivity and/or Ca2+ sensitivity can alter contractile phenotype and electromechanical coupling along the arteriole to capillary continuum. C_LIO_LIThe proposed model can form the basis for detailed multiscale and multicellular models that will investigate cerebral blood flow control in health and in disease. C_LI

Aortic stiffness is both heterogenous and anisotropic. Current non-invasive methods to estimate aortic stiffness are limited to characterizing the aortic tissue as isotropic due to the lack the techniques required to extract multi-axial strain from 3D dynamic images. Vascular deformation mapping (VDM) is a nonrigid image registration technique which has thus far been applied to map aortic growth using longitudinal imaging. In this study, we propose to use VDM to assess 3D aortic deformation by mapping diastolic and systolic images. During image registration process, penalty parameters are employed to fine-tune image alignment and penalize non-physiological deformations. These penalty parameters must be calibrated to ensure that VDM successfully reproduces multi-axial aortic motion patterns in health and disease. In this paper, we developed a calibration pipeline for these parameters using synthetic data. A rotation-free shell model was used to generate physics-based synthetic data on aortic motion incorporating patient-specific geometries, root motion, and blood pressure from a cohort of 14 subjects (healthy, Marfans syndrome and thoracic aortic aneurysm). An error metric was defined to quantify the quality of the VDM results. Furthermore, a k-means clustering technique was used to categorize the subjects into three clusters based on ascending aortic motion. Optimal penalty parameters were identified for each of the three clusters. The results indicated that patient clusters with smaller aortic root motion required larger rigidity penalty values. The calibrated parameters successively reduced errors in 3D displacement and multi-axial stretch compared to un-optimized VDM predictions, enhancing the accuracy of capturing aortic deformation from dynamic images. Among the different aortic regions, the ascending thoracic aorta exhibits the largest error reduction.

Stratified epithelial tissues such as the skin epidermis maintain barrier integrity during development and homeostasis through the coordinated action of cell proliferation, differentiation, delamination, and tissue-scale mechanical forces. During development, the orientation of cell division within the basal layer plays a pivotal role in tissue stratification; however, the mechanical principles linking the orientation of the division plane to these processes across developmental stages remain poorly understood. Here, we expand a recently developed three-dimensional vertex model for stratified epithelia, composed of the basement membrane, basal, and suprabasal layers, to study the mechanical and structural impact of cell divisions with a wider range of orientations. The model integrates developmental stage via specific changes in heterotypic interfacial tensions (arising from actomyosin cortical contractility and adhesion molecules at the basal-suprabasal interface) and tissue stiffness that have been quantified previously in experiments. By systematically varying background mechanical parameters, we investigate how heterotypic tension, division orientation, and tissue fluidity collectively influence the outcome of cell division. Our goal is to uncover the strategies that the embryo may employ to generate stratified phenotypes at different developmental stages, recognizing that these strategies might evolve over time. Although our focus is on the embryonic developmental stages of the epidermis, this framework may also be extended to investigate transformed cells, such as in cancer, to explore how altered division orientation contributes to precancerous or transformed phenotypes.

1Nuclear envelope stretch and rupture are common to cell spreading and migration in a variety of microenvironments, leading to marked changes in nucleocytoplasmic transport. Predicting cell response to different mechanochemical cues that are transmitted to the nucleus remains an open problem in the field of mechanomedicine. We developed a predictive modeling framework to examine how nuclear deformation on substrates with different nanotopographies influences nucleocytoplasmic transport and rearrangement of the nuclear lamina. Using the finite element method, we simulated nuclear compression by the perinuclear actin cap on substrates with arrays of nanopillars, modeling the nuclear envelope as a nonlinear elastic structure and coupling deformations to a biochemical model of lamin remodeling and nucleocytoplasmic transport. These simulations predicted regions of high nuclear envelope stretch adjacent to cell-nanopillar contacts, leading to maximized nuclear envelope tension on small nanopillars spaced by 4-5 microns. We then considered the effects on nuclear transport of YAP and TAZ and found that increased nuclear compression led to YAP/TAZ nuclear localization in agreement with previous experiments. Furthermore, the simulated force load per lamin was maximized on nanopillar substrates with high nuclear stretch. The magnitude of this load was modulated by the rate of actin cap assembly and the overall expression level of lamin A/C - decreasing lamin content in the nuclear envelope led to a higher likelihood of rupture. We validated this prediction in subsequent experiments with lamin-depleted U2OS cells, establishing the central importance of lamin transport and microenvironment nanotopography to nuclear mechanotransduction. 2 SignificanceCell nuclei commonly experience large strains, but existing computational models do not explain the coupling between such deformations and molecular transport. Here, we present a modeling framework that includes the mechanics of nuclear deformations and the reaction-transport of molecules within the cytoplasm, nuclear envelope, and nuclear interior. As a well-controlled setup for comparing experiments and simulations, we consider nuclear indentations exhibited by cells on nanopillar substrates. Our simulations recapitulate measurements of nuclear YAP/TAZ localization from the literature and predict that low-lamin cells experience higher force loads at the nuclear envelope. We validate this prediction experimentally, showing that lamin-depleted cells are more likely to exhibit nuclear rupture. Overall, our framework presents opportunities to predict nuclear mechanoadaptation to different microenvironments.

Pathologic vertebral fractures are a major complication in metastatic spine disease. However, current clinical scores, such as Spinal Instability Neoplastic Score (SINS), show limited predictive capability, particularly within the indeterminate range where most clinical uncertainty lies. This study aimed to develop and evaluate quantitative computed tomography (qCT)-based subject-specific finite element (SSFE) models to predict vertebral strength in presence of different metastatic lesion types. Twelve ex vivo human spine segments, each containing one metastatic (n=12) and one adjacent control vertebra (n=12), were scanned using qCT and calibrated using a calibration phantom. Homogenised nonlinear finite element models were developed with spatially heterogeneous, isotropic, density-dependent material properties and loaded under uniaxial compression corresponding to 1.9% apparent strain. Ultimate failure load, stiffness, and strain distributions were compared between metastatic and control vertebrae. Predicted failure load ranged from 0.2 kN to 6.2 kN (mean. {+/-} standard deviation: 1.8 {+/-} 1.6 kN metastatic; 1.7 {+/-} 1.5 kN control), with no statistically significant difference between groups (p > 0.05). Normalised failure load varied widely, reflecting lesion-specific mechanical heterogeneity. Lytic lesions generally weakened vertebrae, whereas mixed and blastic lesions occasionally enhanced strength, likely due to localised sclerosis or reactive bone formation. High compressive axial strains (greater than 0.019) were frequently concentrated near the endplates, particularly in lytic vertebrae. qCT-derived bone mineral density strongly correlated with failure load (R{superscript 2} = 0.74-0.77). These findings highlight the complexity of metastatic vertebral mechanics and demonstrate that qCT-based SSFE modelling provides a quantitative framework for assessing fracture risk, complementing conventional imaging-based tools.

Cell migration plays a central role in numerous physiological and pathological processes and emerges from the coordinated interplay between intracellular force generation, adhesion dynamics, and mechanical interactions with the environment. A minimal, mechanistically grounded understanding of these processes is required to disentangle the respective contributions of cell-intrinsic and environmental cues. Here, a two-dimensional in silico cell motility model is introduced to describe mesenchymal migration driven by intracellular traction forces generated within actin-rich protrusions anchored to a substrate. The model explicitly accounts for adhesion nucleation, maturation, force buildup and rupture, and relies on a small set of physically interpretable parameters. A systematic mechanical analysis identifies parameter regimes that permit effective cell translocation and delineates conditions leading to stalled or mobile cells. Within motile regimes, the model reproduces a broad spectrum of cell morphologies and migratory behaviours. In particular, cell trajectories exhibit the statistical features of a persistent random walk, with a crossover from ballistic to diffusive motion that arises solely from adhesion dynamics and force balance, without imposing polarization or directional bias. Cell morphology is shown to strongly regulate migration speed, persistence, and pausing behaviour. Altogether, this model provides a minimal reference framework for cell migration on non-deformable substrates and establishes a baseline for future studies of mechanically driven guidance. By construction, it is well suited for extension to deformable fibrous substrates, where cell-induced matrix remodeling and stiffness feedback are expected to bias migration and regulate cell encounters relevant to tissue morphogenesis and anastomosis.

The left ventricle (LV) exhibits torsional deformation during systole, and mechanical relaxation begins during the isovolumic phase. Recent advances in imaging techniques, such as MRI, have revealed that myocardial tissue deformation and sarcomere length changes occur during the isovolumic relaxation phase, even when the chamber volume remains constant. Although such ventricular deformation during the isovolumic phase is considered important for blood ejection and filling efficiency, its mechanistic contribution to contraction and relaxation remains unresolved. In this study, we hypothesized that sarcomere length dynamics during the isovolumic phase affect the isovolumic contraction and relaxation time (IVCT and IVRT) by regulating the contraction force via the force-velocity relationship of ventricular myocytes. To investigate this hypothesis, we focused on experimentally reported differences in the relationship between sarcomere length and LV volume across the endocardial and epicardial layers, as described by Rodriguez et al. We constructed and compared two types of hemodynamic models within the same integrated framework consisting of a circulation model, a LV model, and a myocardial cell contraction model by Negroni-Lascano et al., which differ only in how sarcomere length is determined: a volume-based length model (VL model), in which sarcomere length is uniquely determined by LV volume, and a volume-force-coupled length model (VFL model), in which sarcomere length is determined by the balance between LV volume and contraction force. Simulation results showed that in the VFL model, compared to the VL model, sarcomere length changed during the isovolumic phase, leading to a decrease in contractile force and shortening of IVRT, which may contribute to improved hemodynamic efficiency. These results indicate that sarcomere length dynamics can mechanically regulate force decay during isovolumic relaxation, even under constant left ventricular volume. This study provides a theoretical framework for understanding the contributions of different layers within the LV wall to diastolic function during the isovolumic relaxation phase.

IntroductionArteries, like other soft tissues, exhibit viscoelastic mechanical behavior, meaning their response to stress and strain is time dependent. This implies that the way arteries deform depends not only on the amount of force applied but also on the rate at which the force is applied. This study investigates the effects of different loading rates on the mechanical behavior of human femoropopliteal arteries (FPAs) to understand their rate-dependent characteristics. MethodsHuman FPA specimens were collected from 14 donors, including 7 males and 7 females, aged 45-55 years. A 10x10 mm segment was isolated, mounted onto a biaxial testing device, and subjected to varying loading rates (10 to 50 mN/s). Mechanical responses were recorded, and stress-stretch curves were analyzed. Statistical analyses, including mixed-design ANOVA, assessed the impact of sex and loading rates on tissue stiffness. ResultsResults indicated significant loading-rate dependency, particularly in the circumferential direction. Stretch values decreased with increasing loading rates, more prominently in the circumferential than in the longitudinal direction (p-value<0.01). Statistical analyses revealed no significant interaction between sex and loading rate, though male arteries exhibited slightly higher compliance than female arteries. DiscussionThe findings demonstrate that the mechanical response of FPAs is highly dependent on the loading rate, with more pronounced effects observed in the circumferential direction. At higher loading rates, the human FPAs demonstrated a stiffer response in the circumferential direction. DedicationWe dedicate this work to the memory of our late student, Ali Zolfaghari Sichani, who passed away tragically during his doctoral studies. Ali performed the majority of the experiments and the initial analysis reported in this paper. His passion, dedication, and hard work were the foundation of this research, and he is deeply missed.

Mechanical properties of epithelial tissues play essential roles in morphogenesis and physiological function. In this study, we analytically derived the in-plane bulk modulus, shear modulus, and Poissons ratio of a three-dimensional cell vertex model of epithelial monolayers. We showed that the model can robustly reproduce a near-zero in-plane Poissons ratio, a mechanical feature reported in cultured epithelial tissues. Numerical simulations further confirmed that the theoretically predicted Poissons ratio accurately describes the response of the model under finite, biologically relevant strains. In addition, the model exhibits not only morphological bistability between squamous-like and columnar-like states, but also mechanical bistability characterized by distinct elastic responses. Together, these results provide a minimal three-dimensional framework that links cell-scale mechanical interactions and epithelial morphology to tissue-scale elastic properties.

AO_SCPLOWBSTRACTC_SCPLOWMuscle efficiency decreases with increasing size, largely due to a relative decrease in its mechanical output. Muscle mechanical output depends on its activation, strain, and strain rate and thus varies between different muscles within a limb during locomotion. Distinct muscle coordination patterns are required for efficient cycling, and so we would expect that the coordination patterns for efficient cycling or indeed locomotion would change across animal sizes. We tested whether muscle coordination would change with muscle size using data derived from human cycling: this paradigm allowed for controlled changes in both crank torque and cadence, allowing the multifactorial problem of muscle power output to be decomposed. We used kinematic and pedal data from 12 cyclists undergoing steady pedalling at cadences from 80 to 140 r.p.m. and generated musculoskeletal simulations of their movements. We introduced novel multisegment muscle models in the simulation that incorporated the internal muscle mass and thus accounted for the scaling effects of muscle tissue inertia. We solved the simulations for the muscle activity that was required to minimise the metabolic cost during cycling for each condition. The masses of the muscle models were scaled across five orders of magnitude. The predicted muscle activations were classified by Principal Component analysis to identify whether the coordination of muscle activity was modulated across models with different sized muscles. Analysis of variance revealed significant changes in coordination at the large-scale factors. This study shows how the coordination of muscle activity during locomotion will likely change across a range of body sizes due to the non-linear effects of the inertial mass within the muscle tissues.

Spaceflight-Associated Neuro-ocular Syndrome (SANS) involves complex interactions between intracranial pressure (ICP), intraocular pressure (IOP), and cerebrospinal fluid (CSF) dynamics within the optic nerve subarachnoid space (ONSAS). While existing computational models address specific aspects of these interactions, they lack a comprehensive, system-level representation. To bridge this gap, we present the HEAD (Hemodynamic Eye-brain Associated Dynamics) model. By consistently integrating several previously proposed physiological sub-models, HEAD provides a unified lumped-parameter framework that fully couples cerebrovascular autoregulation, multi-territory ocular hemodynamics, and compartmentalized craniospinal-ONSAS CSF circulation under gravitational loading. This formulation enables the simultaneous analysis of eye-brain-CSF dynamics within a single computational tool. Model predictions were validated against experimental data from supine (0{degrees}) to head-down tilt (HDT, -30{degrees}) postures, accurately reproducing posture-dependent IOP increases and achieving an excellent ICP match against clinical benchmarks at the -6{degrees} HDT standard bed-rest angle. The coupled system predicts bed-specific ocular hemodynamic responses, with retinal blood flow exhibiting the largest relative increase under HDT compared to the ciliary and choroidal circulations. Crucially, explicitly modeling the ONSAS as a distinct compartment reveals a posture-dependent pressure drop of 1.89-3.69 mmHg between the intracranial and perioptic spaces. This compartmentalization yields a translaminar pressure profile that remains positive (8.05-11.83 mmHg) across all simulated conditions but is chronically reduced under sustained HDT. Ultimately, the HEAD model elucidates the physiological mechanisms linking gravitational stress to translaminar mechanics, providing a robust computational foundation to investigate SANS and supply boundary conditions for structural models of the optic nerve head.

Patellar osteochondral allograft (OCA) transplantation is widely used to treat large full-thickness cartilage defects, yet long-term failure and reoperation rates remain high. Although surface congruity and osseous integration are emphasized clinically, cartilage thickness and mechanical compatibility between donor and recipient are not considered. Our previous work suggests that cartilage thickness mismatch can amplify local deformation at the graft boundary, potentially compromising graft longevity. This study investigates how combined mismatches in cartilage thickness and mechanical properties influence the local strain environment at the patellar OCA interface. Simplified two-dimensional axisymmetric finite element models of patellar OCA repair were developed in ABAQUS. Donor-to-recipient cartilage thickness ratios ranging from 0.33 to 3.25 were evaluated together with donor-recipient Youngs modulus mismatches (2.5-7.0 MPa). Cartilage was modeled using homogeneous linear elastic and functionally graded material formulations to account for depth-dependent stiffness. A compressive pressure of 1.0 MPa was applied to represent patellofemoral joint loading, and peak compressive and shear strains were quantified at the graft boundary. Cartilage thickness mismatch produced localized high-strain regions (HSR) of compressive and shear strain at the donor-recipient interface that were absent in thickness-matched constructs. Strain amplification increased with both thickness and mechanical property mismatch. Compressive strain exhibited directional asymmetry, with donor-side-thicker configurations producing greater amplification than recipient-side-thicker configurations. Incorporating depth-dependent cartilage stiffness reduced peak strain magnitudes but did not eliminate mismatch-driven strain amplification. These findings demonstrate that cartilage thickness and mechanical disparity can create HSR at the patellar OCA graft boundary that may predispose grafts to impaired integration and long-term failure.